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Abstract 

We establish a relationship between the zeros of the partition function in the complex mass plane 
and the spectral properties of the Dirac operator in QCD. This relation is derived within the context 
of chiral Random Matrix Theory and applies to QCD when chiral symmetry is spontaneously 
broken. Further, we introduce and examine the concept of normal modes in chiral spectra. Using 
this formalism we study the consequences of a finite Thouless energy for the zeros of the partition 
function. This leads to the demonstration that certain features of the QCD partition function are 
universal. 



1 Introduction 



In recent years the spectral correlators of the Dirac operator in QCD have been the object of intense 
study using both numerical and analytic means. These correlators contain valuable information re- 
garding both the chiral properties of the QCD vacuum and the topological structure of the gauge 
fields. The relation to the chiral properties of the QCD vacuum was established by Banks and Casher 

The eigenvalue density of the QCD Dirac operator at eigenvalue zero is proportional to the chiral 
condensate and is therefore an appropriate order parameter for chiral symmetry. As a complement 
to the Banks-Casher relation, one has the Yang-Lee picture [2] of a phase transition. In the attempt 
to analyze phase transitions in statistical spin models Lee and Yang introduced the concept of the 
zeros of the finite volume partition function in the thermodynamic limit. The volume dependence of 
these zeros allows finite size scaling studies and subsequent identification of universality classes. In the 
case of chiral symmetry, one focuses on the zeros of the partition function in the complex quark mass 
plane. If these zeros pinch the real axis and exhibit a constant density, a discontinuity in the partition 
function arises at the pinch, and chiral symmetry is spontaneously broken. The mass zeros of the 
partition function and the low lying eigenvalues of the Dirac operator thus contain similar information 
about the chiral phase transition. The relation between the two is, however, involved. The partition 
function and its zeros are obtained by averaging over all gauge field configurations. By contrast, the 
eigenvalues of the Dirac operator are given for each gauge field configuration, and only a density of 
eigenvalues is well-defined after averaging. 

The challenge of deriving relations between the zeros of the partition function and the eigenvalues 
of the Dirac operator was first taken up by Leutwyler and Smilga They studied QCD in a Euclidean 
4-dimensional box with side length L subject to the constraint that 

y « L « — . (1) 

Here m.,^ is the pion mass and A is the typical QCD scale. They computed the QCD partition function 
for equal quark masses using the effective chiral Lagrangian and found that quark masses enter only 
in the rescaled combination mL^ S, where S is the chiral condensatein the chiral limit. They further 
observed that the partition function zeros could be thought of as average positions of the eigenvalues. 
While highly suggestive, these results were not completely quantitative. The situation has changed 
dramatically since then. The main break through came with the introduction 4 of random matrix 
concepts in QCD which permits the study of the correlations of the eigenvalues of matrices drawn 
on a general weight constructed to ensure the chiral structure of each eigenvalue spectrum. The 
relation of random matrix theory to QCD in the limit has been established through a number 
of universality studies [3], lattice QCD simulations [HI) [3) a-nd direct calculations using the effective 
chiral Lagrangian |Sj. (For a review of random matrix theory in QCD see 0.) In terms of the spectral 
correlation functions, the universal limit in which QCD and chiral random matrix theory (^RMT) 
coincide is the limit 

iV^oo,A^O,m^O, (2) 

in which the microscopic variables 

C = 2iVA, n = 2Nm (3) 

are kept fixed and A'^ is identified as the dimensionless volume. (Here, A denotes an eigenvalue of 
the Dirac operator and m is the dimensionless quarkmass parameter.) The determination of the 
individual eigenvalue distributions and their most important correlators now permits direct comparison 
of partition function zeros and eigenvalue positions. The suggestion of Leutwyler and Smilga is 
remarkably accurate. The zeros and the average positions of the eigenvalues are intimately connected. 
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Here we shall demonstrate that this relationship can be understood as a fundamental property 
of the chiral ensembles. We show that the zeros are uniquely trapped by the maxima of the joint 
eigenvalue distribution function. This trapping appears on all scales and is thus relevant for any 
finite as well as in the large N limit. To obtain a better understanding of the relation between 
the maxima of the joint eigenvalue distribution function and the zeros of the partition function, we 
introduce and determine the spectral "normal modes" of the chiral unitary ensemble. This provides 
us with a simple tool to describe the fluctuations of the eigenvalues about the maximum of the joint 
eigenvalue distribution. 

As suggested in (0), chiral random matrix theory is not expected to describe all aspects of QCD. 
Only correlations below a certain energy length are expected to be in agreement with chiral random 
matrix theory In solid state physics, this energy is denoted as the Thouless energy. Recently jllj 
it was realized that the effects of a finite Thouless energy can be studied naturally using the language 
of spectral normal modes ^I2j. We thus perform a normal mode analyses of the chiral ensemble to 
formulate and establish certain universal features of the partition function zeros. This argument is 
independent of standard proofs of universality, and its general nature can shed some light on the way 
universality is realized. 

In section |21 we show that for Nf = 1 the zeros of the partition function of xRMT are trapped 
between the maximum positions of the joint distribution function. This result holds for all scales. In 
section El we derive the normal modes of the chiral unitary ensembles and find that they are Chebyshev 
polynomials in the large N limit. We discuss the effects of the Thouless energy in section^ In section 
13 we make the connection to the familiar microscopic spectral density. Our conclusions are contained 
in section ini 



2 Zeros of the partition function in xRMT 
2.1 Chiral random matrix theory 

We start with the partition function of chiral random matrix theory (^RMT) for A'^^ fiavors, which is 
given by |^ El 



Nf 

^ivf/({"i/}) = [ DWll det{D + mf) exp 

f=i 



Ti{W^W) 



(4) 



where (3 denotes the Dyson index and DW is the Haar measure over the Gaussian distributed random 
matrices W . D \s the analogue of the Dirac operator which has the chiral structure 



D 



iW 
iW^ 



(5) 



Here Vl^ is a x M matrix with v = \N — M\ playing the role of the topological charge. Without 
loss of generality we assume v to be positive. The chiral condensate in the chiral limit, S, is related 
to the eigenvalue density of D, p(A), via the Banks-Casher relation ^ 



lim lim lim — ^ 



(6) 



The partition function is invariant under transformations W ^ U'^ WV , where U is a, N x N matrix 
and V a M X M matrix. Following the diagonalization W = WAV, the partition function can be 
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expressed in terms of the eigenvalues of W (with S = 1), 




+00 



+00 



N 

n 

k=l 



N, 



n (Afc + ^ 



-1 

e 



2 -^fc 



A(A 



(7) 



The Vandermonde determinant, A(A^), which is the non-trivial Jacobian of the transformation from 
the matrices to the eigenvalues, has the form 



a(a2) = n (Ai - ) . 

k<l 



(8) 



The partition function Q can now be written as an integral over the joint probability density 
-PAf,^'''(Ai,---,AAr;{m/}) as 



+00 



+00 



N 



n dXkPj 



N,I3 l-^l' 



) A„;{m/}) 



k=l 



(9) 



with 





, TV 


' Nf 




n 


n(Ai + m})Ar^"^e 


\/=i y 


' k=l 





JV/3 
2 



A(A 



(10) 



Unlike real QCD, xRMT has the special feature that the partition function can be expressed in terms 
of the eigenvalues of the Dirac operator. This enables us to derive a number of statements regarding 
the zeros of the partition function. We now focus on the case (3 = 2 — the universality class of QCD 
with 3 colours and quarks in the fundamental representation of the gauge group. (The choice /3 = 1 
corresponds to QCD with two colours in the fundamental representation; /3 = 4 describes QCD with 
any number of flavours and quarks in the adjoint representation of the gauge group.) 

Eq. I|1U() expresses an evident duality between flavor and topology: The joint probability density 
for massless flavours and A'^^ massive flavours depends only on + A'^j,. This relation was proven 
for the QCD partition function independent of xRMT in . 

We now wish to determine the maximum of the joint probability distribution. This will allow us 
to put a tight bound on the zeros of the partition function. The chiral normal modes, to be discussed 
in section 121 describe fluctuations about the maximum of the joint probability distribution. 

2.2 Extremum of the joint probability distribution 

In order to determine the maximum of the joint eigenvalue probability distribution, we consider 
variations of logP^^''^ with respect to the eigenvalues. (We assume the eigenvalues to be ordered with 
Xi < Aj+i.) We introduce the coordinates yi = A? and evaluate the equations 



d log P^^-'" 
dyi 



(11) 



For Nf = and topological sector u this yields 



1 



Nyi 



1 V- 



(12) 
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We now choose to focus on the quenched, i.e. Nf = 0, joint eigenvalue probability distribution. The 
solution to this equation reveals that the maximum of logPjy ' is obtained for 



L-''^ [N\i) = , (13) 

where denotes the generalized Laguerre polynomials. This result follows from the observation that 
Laguerre's differential equation, 

zL%{z)" + (a + 1 - z)L%{zy + NL%{z) = , (14) 

reduces to (|12j) at the zeros of (The proof follows from considerations similar to those made 

in appendix A6 in JH])- For = we can use the fact that 

L-n"\z') = y^^2n(.) (15) 

to see that the density of eigenvalues in the ^ oo limit is precisely that of the usual Gaussian 
ensembles, i.e. a semicircle with support —2 < A < +2. The partition function for /? = 2 is the 
average of a product of fermionic determinants over P^^ 

I Nf TV \ 

<nw}) = (n-rn(^'+-/)) • (16) 
\/=i i=i I 

This can be readily evaluated using orthogonal polynomials, and the result for v = Q agrees with the 
one presented in ^H] 

/=i ^Nf[{-m)}) 



where 



C^^ = det 



</-i(-^-/')l/./^=i,...,^, , (18) 



and Atv^ ^{— m^}^ is the Vandermonde determinant with the negative square of the Nj masses as 
arguments. For the special case of Nf = 1 this yields the result 

zr''"°im) = ^L%' (-iV»^) , (19) 

which is now ready for investigation. The expression (|19|) coincides (up to a constant) with the series 
expansion for the partition function derived in From the expression found there we find that 

Z^'=^'''{m)^m\^\L%\-Nm') . (20) 

2.3 Trapping of the zeros 

The closed forms given above allow us obtain information about the zeros of the partition from the 
spectral correlators. Specifically, we now show that the locations of the partition function zeros in 
xRMT are trapped by the maxima of the joint distribution function. 

The Laguerre polynomials L"{z) are polynomials orthogonal on the interval [0, cxd] with weight 
function w{z) = 2:"exp(— z). They have three properties which are useful for our purpose: 
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• For orthogonal polynomials in general, the zeros of the A^— th order polynomial and the end- 
points of the weight function define intervals. Exactly one zero of the orthogonal polynomial 
of order (A'^ + 1) lies in each of these intervals. 

• For fixed N, the ith zero of the Laguerre polynomial, L^{x) = 0, is a monotonically increasing 
function of a, thus 

^>°- (^') 

• The generalized Laguerre polynomials are related to one another via 

Lt\x) = -'-^^ . (22) 

In the last section we saw that the massless joint distribution function has its maxima when the 
eigenvalues are located at the zeros of L^^^^(iVA^) and that the partition function for Nf = 1 and 
= is proportional to L^{—Nm?). We are interested in relating the zeros of the partition function 
to the position of the eigenvalues at the maximum of the joint distribution function. Since the zeros 
of LP^{—N'm?) follow from those of LP^{N m?) by a rotation from the real to the imaginary axes in m, 
we restrict our attention in the following considerations to Laguerre polynomials of positive argument. 

It follows from (|21j) that the zeros of Z^^ ^''^ ^ are trapped between the corresponding zeros of 
Lj/-^"^ and L^^^^ According to the zeros of l]^"^ lie at the extrema of L^^. The extrema of 
-^AT+i are evidently trapped by the zeros of (and the end points, if necessary). The result is that 

the ith zero of is trapped between the ith and (i + l)st zero of Lj^^'^ and, given the nature of the 
argument, is expected to be closer to the lower value. In other words, the zeros of ^''^ ^{m) are 

trapped by the most probable values for the eigenvalues of P^^^ ^''^ ^ . This result is an exact property 
of the chiral unitary ensembles and is consequently valid on every scale including the microscopic scale. 

The trapping just derived relates the zeros of the A'^y = 1 partition function and the location 
of the maximum in the joint eigenvalue distribution of the quenched {Nf = 0) ensemble. We can 
also relate the zeros of the Nj = 1 partition function to the maximum of its own integration kernel, 

-fjv^ "^''^({Aj}, m). For m — > we have P^-^ ^'^{{^i}^ = Pn^ '^'^"'"^({•^i}) by flavor-topology duality. 
The other limit, m — > oo, decouples one flavor and leaves ^''^ . So, with increasing m, the extrema 
of P^^ ^''^ {m) move smoothly from those of P^^ o,!^+i those of P^^ ^'^ , i.e. from the zeros of 

-^jv ^^'^ ^° ^'^^ zeros of L'^j^^''^. By ((211), they must pass the zeros of the partition function LJ^ for some 
intermediates values of the mass m. The relation between the collective maximum of logP and the 
average eigenvalue positions will be reconsidered in section 121 



3 Normal modes in xRMT 



We have seen that the maximum of the massless joint distribution function is obtained when the 
eigenvalues are located at the zeros of the Laguerre polynomials. In order to study the properties of 
fluctuations about this maximum, it is useful to make a Gaussian approximation to Pjy which 
leads to the form 

log Pi'^=°'" ^ log P^^=°'" + 5\, , (23) 
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where 5Xi is the position of the ith eigenvalue relative to Aj, its value at the collective maximum of 
log Pj^^ ''^ . The matrix C is defined as 

evaluated at the maximum. Concentrating again on the case Nf = 0, we find that the diagonal 
elements of C are 

2u + l ^ (A^ + Af) 
a,^-2iV- — -4gi^. (.5) 



and that the off-diagonal elements are 



We now consider the eigenvalue equation for the real symmetric matrix C: 

N 



Y^C,,<pf = u:,<^f\ (27) 
1=1 



The eigenvectors, (f)^ > , are the (normalized) normal modes of the xRMT spectrum. They describe the 
statistically independent correlated fluctuations of the eigenvalues of the random matrix about their 
most probable values. The normal modes provide an alternate description of the eigenvalues of any 
given random matrix since 

AT N 

5\i = cuc^f^ with ^?^f'^ = ^kk' . (28) 

k=l i=l 

We can locate the eigenvalues by specifying either the dXi or the amplitudes Ck as convenient. The 
eigenvalues, w^, provide a measure of the magnitude of these fluctuations. 

The derivation of these eigenvalues and eigenvectors can be performed as in jE]. The resulting 
eigenvalues are 

ujk = -AkN. (29) 

As in "12_' , we flnd a linear dispersion relation valid for all k and N. The linearity of ()29|) is a reflection 
of the well-known rigidity of random matrix spectra. Furthermore, this result is independent of v. 
Just as in the Gaussian case, the eigenvectors are found to be Chebyshev polynomials in the large N 
limit (i.e., with corrections of order 1/N): 

#^ = V^^^^-i(y) • (30) 

The normalization of the eigenvectors is 

^ dxp{x)<P^''Hx)^^'Hx) = 6ki , (31) 

where p{x) is again the semicircle 

N 



p{x) = — \/4-x2 . (32) 



TT 
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Note that only odd normal Chebyshev polynomials appear. This is a consequence of chiral symmetry, 
which ensures that all non-zero eigenvalues come in pairs (—A, A). Equation (|23|) can now be written 
as 

'piV/=0,i/\ ^ N 

' -^El^fcl'^'^- (33) 



log 



pNf=0,u 



k=l 



Following (|28|) . the coefficients Ck are constructed as 



N 

Ck = Yl ^f"* ' 
i=l 



(34) 



and are statistically independent. It is obvious from that that the mean square amplitude for the 
A;th normal mode is 



(4) 



1 



(35) 



We can use the normal modes to construct a Gaussian approximation to the partition function as 



■'N 



(m) 



N 



Y[ dck exp 



k=l 



N 



k=l 



(m2 + (Xk + SXk)^) 



(36) 



with the Afc given by ((T^ and 5Xk given by (^5]) . The Gaussian approximation also permits a 
simple approximate calculation of the number variance jl2j . This calculation reveals that the familiar 
logarithmic behaviour of the number variance (i.e., the "spectral rigidity" of the random matrix 
ensembles) is a direct consequence of the linearity of the dispersion relation H29|) for all k. 

We have chosen to consider the normal modes for the case Nf = and z/ = 0. This choice is some 
what arbitrary; it would be equally sensible to start with the Nf = and = 1 normal modes. We 
will employ this Gaussian approximation below to consider the sensitivity of the partition function 
zeros to the effects of a Thouless energy. Since the resulting shifts are small, this arbitrariness will be 
of no consequence. 



4 Effects of a Thouless energy 

Normal modes describe the correlated fluctuations of eigenvalues about their most probable values. 
As we have seen ()29() . the normal modes for ^RMT obey a linear dispersion relation. By contrast, 
uncorrelated eigenvalues obey a quadratic dispersion relation, and the mean square amplitude of the 
lowest mode with A; = 1 is larger by a factor of N/k In QCD, it is expected that spectral 

correlations in a sufficiently small energy domain will follow the predictions of xRMT. On larger 
energy scales, spectral correlations die out. The characteristic energy which divides these regions is 
the Thouless energy, usually denoted by Ec- In 4-dimensional QCD the Thouless energy is estimated 
to be [im 

EJD ~ y/N , (37) 

where D is the mean level spacing. This behaviour has been verified in lattice studies |18j . The 
connection between the Thouless energy and the normal modes of the eigenvalue spectrum has been 
investigated in for the case of sparse matrices. There it was found that "almost all" normal 
modes obey the linear dispersion relation discussed above with remarkable accuracy. Reflecting the 
presence of a Thouless energy, a small fraction (i.e., approximately 1/y/N) of long wave length modes 
are strongly enhanced. The mean square amplitude of the longest wave length mode with k = 1 
approaches the value appropriate for uncorrelated eigenvalues. Since such enhancement can be readily 
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incorporated in our Gaussian approximation to the partition function (|36j) . normal modes provide 
us with a natural and convenient tool to study the effects of a Thouless energy on the zeros of the 
partition function. So far we have seen that the most probable eigenvalues are located at the zeros of 

— 1/2 

the Laguerre polynomial and that the zeros of the partition function are given by the zeros of 

L^. The Gaussian approximation allows us to see how this result is modified by the presence of 
a Thouless energy. To mimic the effects of the Thouless energy, we shall enhance the long wavelength 
modes in the partition function. 

Our aim is to demonstrate that the zeros of the partition function in the microscopic region remain 
virtually unaffected even if the enhancement of the soft modes is substantial. Since we are concerned 
only with the microscopic zeros, every long wavelength mode contributes to a "breathing" of the 
spectrum. In order to investigate the influence of the soft modes, it is sufficient to evaluate the 
strength (i.e., mean square amplitude) of this effective breathing. 



For concreteness we start with yN longest wavelengths modes with fluctuations as given by the 
Gaussian approximation. Let Aj denote the values at the maximum of ''^ . The fluctuations 



introduced in the Gaussian approximation by the \/N longest wavelengths modes are 

~Xi^\i\l + Y. r?(^') Ck\=hs . (38) 

k=l 



From pUj) we know that the normal modes are \J2/N U2k~i i^f^)- This gives for the linear coefficient 

r/(') = /|M-l)'+' • (^^^ 

Additionally we know from (|26|) that the normal modes are linearly independent. With this information 
the variance of s becomes in leading order of 

k=l k=l 

Here we used the facts that the linear terms vanish and that (c^) = with (|29|) for the quadratic 

terms. The result is simply a (D{l/N) correction. 

In order to study the effects of a Thouless energy, we now enhance the mean square amplitudes 
of these y/N soft modes by a factor of N/k"^. This factor provides a smooth interpolation from the 
behaviour of uncorrelated soft modes (for A; = 1) to that of xRMT (for k = \/iV)- This interpolation 
is completely consistent with the results of We now find that 

/,2. ..2 1 N log(jV) 

k=l 

The result still shows strong suppression in N . The decision to single out ^/N soft modes for enhance- 
ment is not essential. We could declare any fraction of the long wave length modes which vanishes 
in the N ^ oo limit as "soft" . A similar interpolation between the limits of uncorrelated and ^RMT 
modes will always lead to a value of {s"^) — 1 which vanishes as ^ oo. In short, the effects of a 
Thouless energy are expected to have a negligible effect on (A?) over the entire microscopic spectrum. 

The question is now how we can evaluate the effect on the zeros from the enhanced long wave 
length modes. To this end we introduce the distribution function of the fluctuations 

PNis) = ^s^e-^ (42) 



9 



and fix 7 and cr by the value of (s^) found above. M is the normalization. In the RMT case we have 
7 = N/2 — 1 and cr^ = 2(A^ + 1)/N'^, while in the case where the long wavelength modes are enhanced 
we have 7 = iV/log(iV) - 1 and a'^ = log(iV)(l + log(iV)/(27V))/iV. 

The effect of the first ^/N normal modes is different for large and small eigenvalues. Whereas for 
the smallest ^/N eigenvalues it amounts to a breathing it means incoherent fluctuations for the larger 
ones. We now evaluate the effect of this breathing on the microscopic zeros. Recall that the partition 
function for Nf = 1 and = is the average of the fermionic determinant with respect to the joint 



probability distribution, and that H 
quantity = 2 iV m is fixed we have 



i(m2 - A2) = L~^'^{Nm'^). In the microscopic limit where the 



1 



]S[\22N 



H2N 




(43) 



To investigate the correction to the microscopic partition function zeros under the influence of the 
y/N longest wavelength normal modes we thus have to consider the following integral 



1 f°° ( \ 



cos 



where the normalization factor /C is 



/C 



ds PAr(s) S 



2VN 



(44) 



(45) 



The factor s^^^ comes from the rescaling in the fermionic determinant. For the evaluation of (|44|) we 
make use of the series expansion of the cosine 



X 

cos I — 
s 



cos(x) + X sin(x) (s — 1) + 0{s — 1) 



This yields 



Z{n) ~ cos(/u) + fism{fi) 



V2(jr((7 + 2ViV + 2)/2 



r((7 + 2^+l)/2 



+ 



(46) 



(47) 



Since 7 \/iV, we can use that [TO] 

T{z + g) 
T{z + b) 



a~b 



1 + 



(a-6)(a + 6- 1) 



2z 



+ 



(48) 



and finally find that in the case where we consider v N RMT long wavelength modes the correction 
term to cos(/u) is of order l/\/iV, and in the case of \fN enhanced wavelength modes it is of order 
\ogN/y/N . In both cases the zeros of the partition function are unaffected in the microscopic limit. 

These results suggest a new kind of "universality" of the microscopic partition function. As 
a specific example this universality shows that the microscopic zeros are unaffected even if the \/iV 
longest wavelength normal modes are enhanced in such a way that they interpolate between Poissonian 
and RMT statistics. 



5 The microscopic limit 

So far we have been discussing the joint distribution function, the zeros of the partition function, 
and the N positions that specify the collective maximum. Here we link this to the more familiar 
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microscopic-eigenvalue density. For a finite the eigenvalue density is found from the joint distribu- 
tion function as 



p/' (A,{m/})=/ .../ dX2...dXNPN'iX,X2,h,---,^N;{mf}). (49) 

J —oo J —oo 

The double-microscopic spectral density is then defined as [U 1201 

and similarly for all other spectral correlations. The functional form of the microscopic eigenvalue 
density has been derived in [2j| and for Nf = 1 and topological charge v it reads 



Ps' {C:p) = ^[Ju-i{C) +Ju-2{C)UC))-\C\ (c^ + ^2)j^_,(^) • (^1) 

The partition function in the microscopic limit for one flavor is proportional to luilj)- The extrema 
of (jKT|) are given by 

dC 

and are obviously functions of /x. If the derivative in (|52|) is evaluated at = /i then the solutions 

of (|52|) corresponding to the maxima of ph^^ ^''^^ coincide with the zeros of the Nj = 1 microscopic 
partition function, Z^f^^''^{m), for an imaginary argument. To conclude: The peaks of the microscopic 
eigenvalue density are in one to one correspondence with the zeros of the partition function and by the 
trapping proven in section ITITI also to the positions which maximize the joint probability distribution. 



6 Conclusions and outlook 

In the present paper we have established an intimate relationship between zeros of the partition 
function and the spectral properties of the Dirac operator. The relation is derived within chiral 
random matrix theory and applies to QCD Dirac spectra and partition function zeros near the origin 
(in the microscopic regime). Through the introduction of spectral normal modes we have tested the 
validity of the relationship when a finite Thouless energy is introduced. The observed independence 
of the microscopic zeros compliments the existing universality studies. The present study treats one 
flavour. For more flavours with degenerate masses, it is known that the zeros of the microscopic 
partition function are not confined to the imaginary axis [22]. While this makes the relation between 
the eigenvalues and the zeros somewhat less direct, there is no reason to expect that the normal mode 
analyis should not apply for any number of flavours. 

We remark that the normal mode analysis is generically applicable and not a special feature of 
random matrix theory. In particular, the normal mode analysis lends itself to a study of the spectral 
properties of the Dirac operator in lattice QCD. Such a study would be truly interesting in that it 
would shed new light on the role of random matrix like correlations in lattice gauge theories. Almost 
all of the normal modes are expected to be given by random matrix theory while only the very long 
wavelength modes are determined by the detailed dynamics. 

In a broader perspective this study may also be seen as the first step towards the establishment of 
a one-to-one correspondence between the zeros and the most likely eigenvalue positions whenever the 
short range spectral correlations are random matrix like. If such a general relation were established 
then it could be used to argue that the critical exponents for the Yang-Lee edge in QCD must coincide 
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with the ones for the gap in the spectral density of the Dirac operator. Such relationship, if true, would 
allow for substantial simplifications when trying to determine critical exponents by lattice techniques. 
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